Constraining the dark energy and smoothness-parameter with supernovae 
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The presence of inhomogeneities modifies the cosmic distances through the gravitational lensing 
effect, and, indirectly, must affect the main cosmological tests. Assuming that the dark energy is a 
smooth component, the simplest way to account for the influence of clustering is to suppose that the 
average evolution of the expanding Universe is governed by the total matter-energy density whereas 
the focusing of light is only affected by a fraction of the total matter density quantified by the a Dyer- 
Roeder parameter. By using two different samples of SNe type fa data, the Q m and a parameters are 
constrained by applying the Zeldovich-Kantowski-Dyer-Roeder (ZKDR) luminosity distance redshift 
relation for a flat (ACDM) model. A x 2 -analysis using the If 5 SNe la data of Astier et al. sample 
(2006) constrains the density parameter to be Q m — 0.261q 'Jy (2<t) while the a parameter is weakly 
limited (all the values € [0, f] are allowed even at la). However, a similar analysis based the f82 SNe 
la data of Riess et al. (2007) constrains the pair of parameters to be f2 m = 0.33^0 07 and a > 0.42 
(2a). Basically, this occurs because the Riess et al. sample extends to appreciably higher redshifts. 
As a general result, even considering the existence of inhomogeneities as described by the smoothness 
a parameter, the Einstein-de Sitter model is ruled out by the two samples with a high degree of 
statistical confidence (f f .5a and 9.9a, respectively). The inhomogeneous Hubble-Sandage diagram 
discussed here highlight the necessity of the dark energy, and a transition deceleration/accelerating 
phase at z ~ 0.5 is also required. 

PACS numbers: Dark energy, cosmic distance, supernovas, inhomogeneities 



I. INTRODUCTION 

The Hubble-Sandage diagram for Type la Supernovae 
(hereafter SNela). as measured by the Supernova Cos- 
mology Project [1| and the High-z Supernova Search 
Team[2[, provided the first evidence that the present Uni- 
verse is undergoing a phase of accelerating expansion 
driven by an exotic component with negative pressure 
(in addition to the cold dark matter), usually called dark 
energy. 

The idea of a dark energy-dominated universe is a di- 
rect consequence of a convergence of independent ob- 
servational results, and constitutes one of the greatest 
challenges for our current understanding of fundamental 
physics [3(. Among a number of possibilities to describe 
this dark energy component, the simplest and most the- 
oretically appealing way is by means of a cosmological 
constant A, which acts on the Einstein field equations 
as an isotropic and homogeneous source with a constant 
equation of state, w = p/p = — 1. 

Although cosmological scenarios with a A term might 
explain most of the current astronomical observations, 
from the theoretical viewpoint they are plagued with at 
least a fundamental problem, namely, it is really difficult 
to reconcile the small value of the vacuum energy den- 
sity required by observations (~ 10~ 10 erg/cm 3 ) with es- 
timates from quantum field theories ranging from 50-120 
orders of magnitude larger This problem sometimes 
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called the cosmological constant problem (PCC) has in- 
spired many authors to propose decaying A models Q 
and other alternative approaches for describing dark 
energy Q. Nevertheless, the present cosmic concordance 
model (CCM) which is supported by all the existing ob- 
servations is a flat ACDM cosmology with a matter frac- 
tion of f2 m ~ 0.26 and a vacuum energy contribution of 

n A ~ 0.740, Si, E3. 

On the other hand, the real Universe is not perfectly 
homogeneous, with light beams experiencing mass inho- 
mogeneities along their way thereby producing many ob- 
servable phenomena. For instance, light lines traversing 
in the universe are attracted and refracted by the gravita- 
tional force of the galaxies on their path, which bring us 
the signal of lensing, one of which is the multiple images 
of a single far galaxy [Til [T^. Nowadays, gravitationally 
lensed quasars and radio sources offer important probes 
of cosmology and the structure of galaxies. The optical 
depth for lensing depends on the cosmological volume el- 
ement out to moderately high redshift. In this way, lens 
statistics can in principle provide valuable constraints on 
the cosmological constant or, more generally, on the dark 
energy density and its equation of state [Tj, Q EH ■ 

In this context, one of the most important issues in 
the modern cosmology is to quantify from the present 
observations the influence of such inhomogeneities on the 
evolution of the Universe. An interesting possibility to 
account for such effects is to introduce the smoothness 
parameter a which represents the magnification effects 
experienced by the light beam. When a = 1 (filled 
beam), the FRW case is fully recovered; a < 1 stands for 
a defocusing effect; a = represents a totally clumped 
universe (empty beam). The distance relation that 
takes the mass inhomogeneities into account is usually 
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named Dyer-Roeder distance [l6j. although its theoretical 
necessity had been previously studied by ZeldovichJT?} 
and Kantowski[l8j]. In this way, we label it here as 
Zeldovich-Kantowski-Dyer-Roeder (ZKDR) distance for- 
mula (for an overview on cosmic distances taking into 
account the presence of inhomogeneities see the paper 
by Kantowski[l9(). 

Several studies involving the ZKDR distances in dark 
energy models have been published in the last few years. 
Useful analytical expressions for ACDM models have 
been derived by Kantowski et al. [2(1, HH and Demianski 
et a/.(22j. Working in the empty beam approximation 
(a = 0), Sereno et al.^j^ investigated some effects of the 
ZKDR distance for a general background. By assuming 
that both dominant components may be clustered they 
also discussed the critical redhift, i.e., the value of z for 
which dji{z) is a maximum (or Q(z) minimum), and com- 
pared to the homogeneous background results as given by 
Lima and Alcaniz [241 ] , and further discussed by Lewis and 
Ibata[25j. Demianski and coworkers derived an approx- 
imate solution for a clumped concordance model valid 
on the interval < z < 10. Additional studies on this 
subject is related to time delays(2f| [26j], gravitational 
lensing[13, HI], and even accelerated models driven by 
particle creation have been investigated]!^]. 

In a previous paper[3(|, we have applied the ZKDR 
equation in the framework of phantom cosmology in or- 
der to determine cosmological constraints from a sam- 
ple of milliarcsecond compact radio sources. By assum- 
ing a Gaussian prior on the matter density parameter, 
i.e., O m = 0.3 ±0.1, the best fit model for a phan- 
tom cosmology with u> = —1.2 occurs at Q m = 0.29 
and a — 0.9 when we marginalize over the characteristic 
size of the compact radio sources. Such results suggested 
that the ZKDR distance can give important corrections 
to the so-called background tests of dark energy. In this 
article, the pair of cosmic parameters, Q m = 1 — J7a 
and a, are constrained from Supernovae observations by 
applying the ZKDR distance-redshift relation for a flat 
(ACDM) model. As we shall see, the a parameter is 
not well constrained by the 115 SNe observed by Astier 
et al. (2006). However, the 182 SNe Type la sample 
of Riess et al. (2007), constrains the pair of parame- 
ters to be Q m = 0.33t°;°7 and a > 0.42 (2a). As a 
general result, even considering the existence of inhomo- 
geneities described by the a parameter, the Einstein-de 
Sitter model is ruled out by the two samples with a high 
degree of statistical confidence (11.5a and 9.9c, respec- 
tively) . 

The paper is organized as follows. In section II, we 
present the basic equations and the distance description 
taking into account the inhomogeneities as described by 
the ZKDR equation. In section III, we determine the 
constraints on the cosmic parameters from the two Su- 
pernovae samples. Finally, we summarize the main con- 
clusions in section IV. 



II. ZKDR EQUATION FOR LUMINOSITY 
DISTANCE 

In a clumpy universe model, the local geometry is in- 
homogeneous, but its global aspect can be described by 
the FRW type geometry (c = 1) 

ds 2 = dt 2 - R 2 (t)(dr 2 +r 2 dn 2 ), (1) 

where R(t) is the scale factor and dVl 2 denotes the metric 
in the 2-sphere. 

As it is widely known, the idea of clumpy universe is 
still a ill-defined notion since we do not have a clear math- 
ematical recipe to separate the global properties from the 
local inhomogeneous aspects of the Universe. After Dyer 
& Roeder[l6j], it is usual to introduce a phenomenologi- 
cal parameter, a = 1 — < ^ ci > , called the "smoothness" 
parameter. Such a parameter quantifies the portion of 
matter in clumps (p c i) relative to the amount of back- 
ground matter which is uniformly distributed (p m ). In 
general, due to the structure formation process, it should 
be dependent of the redshift, as well as, on the direction 
along the line of sight (see, for instance, [l4l.[30| and Refs. 
there in). However, in the majority of the works a is as- 
sumed to be a constant parameter. From a mathemati- 
cal viewpoint the treatment is based on the optical-scalar 
equation for light pro pag ation in the so-called geometric 
optics approximation pll [3l| 

VA" + -R pv k»k v VA = 0, (2) 

where a prime denotes differentiation with respect to the 
affine parameter A, A is the cross-sectional area of the 
light beam, i? M „ the Ricci tensor, and the photon 
four-momentum. In this form, it is implicit that the in- 
fluence of the Weyl tensor (shear) can be neglected. This 
means that the light rays are propagating far from the 
mass inhomogeneities so that the large-scale homogeneity 
implies that their shear contribution are canceled. The 
proportionality factor between the cross-sectional length 
A? and the angular distance dA can be defined to be 
constant. Actually, the above optical scalar equation 
is usually written in terms of the dimcnsionless angu- 
lar diameter distance Da = HgdA- Further, by recall- 
ing the existence of a simple relation between the lumi- 
nosity distance, and the angular-diameter distance (from 
Etherington principle [32|. dr, = (1 + z) 2 cU), it is easy to 
show that the ZKDR (dimcnsionless) luminosity distance 
for ACDM cosmology satisfies the following differential 
equation [H HH, HI 

(l + z fF^-(l + z )g^±+HD L = 0, (3) 

which satisfies the boundary conditions: 
D L (0) - 0, 

(4) 

dD L I _ -, 

STlo - 1- 
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FIG. 1: The a-effect on the residual magnitudes. In Fig. a we show the 115 Supernovae data from Astier et oZ.0], and the 
predictions of the ZKDR luminosity distance for several values of a relative to an empty model (fi m = 0, Qa = and a = 1). 
In Fig. b we show the same graph but now for the 182 SNe type la from Riess et al. sample [§[. For comparison, in both 
panels we see (black curves) the prediction of the cosmic concordance model (Q m = 0.26, = 0.74, a = 1). 



where T ', Q end Ti are functions of the cosmological pa- 
rameters, expressed in terms of the redshift by: 

t = n rn + (i-a m )(i + z)- 3 , 

Q = ^ + 2(l-fl m )(l + z)- 3 , 

n = (^ Y ^n m + 2(i-a m )(i + z )- 3 , (5) 

where as remarked before, the a parameter appearing in 
the TL expression (here assumed to be a constant) quan- 
tifies the clustered fraction of the pressureless matter. 

III. SAMPLES AND RESULTS 

The standard FRW models contain only homoge- 
neously and isotropically distributed perfect fluid gravity 
sources, and the present CCM is assumed to represent 
both the "large scale" geometry of the universe and the 
matter content. However, the Universe appears homoge- 
neous only in a statistical sense, when one is describing 
the largest scales. Therefore, although making very use- 
ful predictions, our cosmological models are somewhat 
inadequate at small and moderate scales. This means 
that relations like ^(Hq, Q m , A; z), the distance modulus 
for a standard candle, commonly assumed to be valid 
on average can be incorrect even for observations includ- 
ing SNe la. In particular if the underlying mass density 
approximately follows luminous matter (i.e., associated 
with bounded galaxies), the effects of inhomogeneities on 
relations like ^(Hq, Q rn , A, z) must be taken into account. 



In Fig. 1, we display the effects of the inhomogeneities 
in the reduced Hubble-Sandage Diagram for the Astier 
et al. (2006) and Riees et al. (2007) samples for some 
selected values of the smoothness parameter. The plots 
correspond to several values of Q m and a as indicated in 
the panels. The difference between the data and mod- 
els from an empty universe case (OCDM) prediction is 
also displayed there. For the sake of comparison, we also 
show the Einstein-de Sitter (E-dS) model, i.e. il rn = 1 
and a = 1, as well as the present cosmic concordance 
(Q m = 0.26, fl A = 0.74, a = 1). Note that the a pa- 
rameter contributes in the right direction i.e., the SNe 
type la become dimmer when it increases on the allowed 
range. In what follows, a y 2 minimization will be ap- 
plied for the two sets of SNe data with the parameters 
f2 m and a spanning the interval [0,1] in steps of 0.01, for 
all numerical computations. 



A. Astier et al. Sample (2006) 

Let us now discuss the bounds arising from SNe la 
observations on the pair of parameters (Q m ,a) defining 
the ZKDR luminosity distance. 

The current data from Supernova Legacy Survey 
(SNLS) collaboration correspond to the first year results 
of its planned five year survey. The total sample includes 
71 high-z SNe la in the redshift range 0.2 < z < 1 plus 44 
low- z SNe la as published by Astier et aL0- Although in 
a better agreement with WMAP 3-years results Q than 
the gold sample @| (for a more detailed discussion see 
e.g., Jassal et aL[33|), the most distant SNe la of these 
115 events has redshift smaller than unity. 

Following standard lines, the maximum likelihood es- 
timator, CsNia oc exp [— x%Nia( z > p)/2] 7 is determined 
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FIG. 2: a) The f2 m — a plane for flat ACDM models obtained from 115 SNe la data Astier et aZ.Q. Note that the a parameter 
is not well constrained by the data, b) The Likelihood for the a smoothness parameter. We see that even at la the smoothness 
parameter is poorly restricted (all its admissible values are allowed), c) Probability of the matter density parameter. We see 
that 0.19 < f2 m < 0.43 with 2a confidence level. 



by a x 2 statistics 

2 / i \ \ - (m(z»;p) ~ ^o,») 2 f(l s 

i Mo,< int 

where p = (H , a, O m ) is the complete set of parameters 
that we want to fit, a a „ , , cr !tlt are, respectively, the errors 
associated with the observational techniques in determin- 
ing the distance moduli (includes a peculiar contribution) 
and the intrinsic dispersion of SNe la. The corresponding 
errors are reported in Astier et al. paper Q. 

Marginalizing our likelihood function over the nuisance 
parameter, Hq, we obtain the likelihood function for the 
&>m — ol plane. In order to determine the cosmological 
parameters (Cl m ,a), a \ 2 minimization for the range of 
[0,1] in steps of 0.01 has been applied. The 68.3%, 90.0% 
and 95.4% confidence levels are defined by the conven- 
tional two-parameters x 2 levels 2.30, 4.61 and 6.17, re- 
spectively. It is very important to note that we do not 
consider any prior in fl m , as usually required by the SNe 
la test. The basic results are shown in Figs. 2a, 2b and 
2c. From fig. 2a we see that all range for a is accepted, 
while a S7 m ss 0.3 ±0.1 is obtained. In Fig. 3b we see the 
likelihood for the smoothness parameter. The best fit ad- 
justment occurs for values of a = 1.0 and fi m = 0.26 with 
Xmin = H3.3 and v — 113 degrees of freedom (x 2 l v = 1), 
thereby showing that the model provides a good fit to 
these data. It is also interesting that for any a value, we 
also find no evidence for a high Q m parameter as required 
by a flat Einstein-de Sitter universe (ft a = 0). Actually, 
the Einstein-de Sitter scenario has a very small statistical 
significance x 2 = 244.9 (11. 5c outside). However, since 
the Astier et al. data are not restrictive for the a pa- 
rameter, let us now consider the enlarged SNela sample 
observed by the High-z Supernovae Search Team Q. 



B. Riess et al. Sample (2007) 

The so-called gold sample from the HZS team Q is a 
selection of 182 SNe la events distributed over the red- 
shift interval 0.01 < z < 1.755, and constitutes the com- 
pilation of the best observations made so far by them 
and by the Supernova Cosmology Project events observed 
by Hubble Space Telescope (HST). As before, constrains 
on the cosmological parameters (fl m ,a), are determined 
from a % 2 minimization within the range of [0,1] spanned 
by such parameters. 

In fig. 3a one can see that 0.42 < a < 1.0 and 
0.25 < O m < 0.44 with 90% of statistical confidence. The 
best fit adjustment occurs for values of fi m = 0.33 and 
a = 1 with Xmin = 158.6 and v — 180 degrees of freedom 
the reduced (x 2 '' /v ~ 0.9). Therefore, the model provides 
a very good fit to the Riess et al. sample. In Fig. 3b we 
see the likelihood for the smoothness parameter. As pre- 
viously remarked, the Riess et al. data set is much more 
restrictive for the smoothness parameter than Astier et 
al. sample. Within 2a, the allowed range for the a falls 
on the interval 0.42 < a < 1.0 (cf fig 2a). In fig 3c we 
show the probability for the density matter parameter. 
In this case a small region is permitted 0.25 < fl m < 0.44 
with of the confidence level (2a) . In the analysis for the 
Einstein-de Sitter universe (fi ro = 1.0, f^A = 0.0), the 
X 2 = 255.8 is too bad (9.9er c.l. outside for 1 degree of 
freedom), and guarantee us to exclude this model with 
high confidence. 



IV. COMMENTS AND CONCLUSIONS 

Cosmology is in an exciting period. A considerable 
set of rather sophisticated experiments, until a few years 
ago regarded as futuristic, have now been completed with 
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FIG. 3: a) Confidence contours on the (Q m ,a) plane for flat ACDM models as inferred from 182 SNe la measurement by 
Riess et al. 0. b) The Likelihood function for the a smoothness parameter. We see that at 2a the smoothness parameter is 
restricted on the interval (0.42 < a < 1.0). c) Probability of the matter density parameter. In this case a comparatively small 
region is permitted 0.25 < Q m < 0.44 with (2a) confidence level. 



spectacular success. The results of the first observations 
almost one decade ago have been confirmed what was 
long surmised, namely, that most of the matter is non- 
baryonic and that we live in an accelerating expanding 
Universe dominated by dark energy. 

In this article, we have discussed the influence of in- 
homogeneities on the expansion rate of the universe. In 
particular, if the smoothness a parameter could be con- 
strained through a statistical analysis involving two large 
sets of SNe la data. As we have seen, in the case of the 
Astier et al. sample, the entire interval of a is allowed 
while a fl rn w 0.3 ± 0.1 is obtained. Within the existing 
uncertainties, these results are consistent with the con- 
straints obtained from angular diameter of compact radio 
sources with basis on the Gurvits et al. data [3(| HH). 
Therefore, although in close agreement with rather differ- 
ent analysis, this SNe data set is uncapable to constrain 
the smoothness parameter. Actually, at this moment, 
the sample of Riess et al. provides a more stringent con- 
strain with the allowed range for a falling on the inter- 
val 0.42 < a < 1.0 (2a). Basically, this occurs because 
the Riess et al. sample extends to appreciably higher 
redshifts. In general, both analysis suggest that a large 
range for a is permitted, and that the Einstein-de Sitter 
model is strongly excluded (11.5er and 9.9er, respectively). 
As we have seen, the necessity of the dark energy and a 
transition from deceleration to an accelerating phase is 
maintained even when one takes into account the clus- 
tering phenomenon. However, at the level of the SNe 
la observations discussed here, these results suggest that 
the dumpiness of matter distribution can mimic at least 



a small fraction of the dark energy component. 

Finally, we would like to stress that measurements 
from SNe la combined with the ZKDR inhomogeneous 
approach adopted here may provide an independent and 
more rigorous cosmological test for the cosmic concor- 
dance model in the near future. In this concern, it 
should be very important to investigate whether the a 
parameter can be constrained using independent obser- 
vations, among them: the cosmic microwave background 
anisotropies, the physics of galaxies clusters, Sunyaev- 
Zeldovich effect, time delay and statistical of gravita- 
tional lensing. Some studies along these lines will be 
presented in a forthcoming communication. The present 
results based only on the Hubble-Sandage diagram show 
that the Riess et al. sample more restrictive than the 
Astier et al. sample, thereby reinforcing the interest to 
observe more supernova events at higher redshifts. 
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